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1. Introduction 

The motivation of this work comes from the problem of numerical modeling of acoustic scattering in 
tubular waveguides geometrically perturbed up to infinity (e.g. bent or stretched) and filled with homoge- 
neous or inhomogeneous media. In order to obtain a good approximation of scattered waves by numerical 
solutions of a problem with finite computational domain, waveguides should be truncated without creating 
excessive reflections from the artificial boundary of truncation. The idea is to place in front of the boundary 
of truncation a layer strongly absorbing the scattered waves. Due to the strong attenuation, the homoge- 
neous Dirichlet boundary condition is a suitable boundary condition on the boundary of truncation. This 
truncation scheme supplemented with very special construction of the layer is widely known as the Perfectly 
Matched Layer (PML) method, originally introduced in pQ. The method is in common use for numerical 
analysis of a wide class of problems. For some of them, stability and convergence of the method have been 
proved mathematically; e.g. [2] [7J [TTJ [14] . In the present paper we develop the PML method for Neu- 
mann Laplacians on manifolds with axial analytic quasicylindrical ends and prove stability and exponential 
convergence of the method. Neumann Laplacians model the scattering problem described in the beginning 
of introduction, see, e.g., J3J. 

As is known, construction of PMLs is closely related to complex scaling. Complex scaling involves complex 
dilation of variables and has a long tradition in mathematical physics and numerical analysis [4j [16] . In this 
paper we construct PMLs in a different way. Instead of complex dilation of variables we replace the metric on 
a part of the manifold by a non-degenerate complex symmetric tensor field. This approach is close in spirit 
to |14] and can be understood as a deformation of the Remanian geometry by means of the complex scaling. 
As a result of this deformation all formulas for the quadratic form and for coordinate representations of the 
Neumann Laplacian turn into the corresponding formulas for the quadratic form and the non-selfadjoint 
operator describing infinite PMLs. This essentially simplifies construction and tractability of the formulas. 
Let us stress here that due to variation of the metric along quasicylindrical ends not only the Laplace- 
Beltrami operator but also the Neumann boundary condition should be changed in PMLs. This new effect 
leads to significant difficulties in analysis of the PML method. 

Relying on ideas of the Aguilar-Balslev-Combes-Simon theory of resonances [HUB] we establish a limiting 
absorption principle. As is typically the case, scattered solutions satisfying the limiting absorption principle 
locally coincide with solutions to the problem with infinite PMLs. The latter solutions are of some exponential 




Figure 1. Representation M = M c U II. 



decay at infinity. These results are mainly based on localization of the essential spectrum of non-selfadjoint 
operators corresponding to the problem with infinite PMLs. Thanks to the exponential decay of solutions 
in PMLs we can establish stability and exponential convergence of the PML method by using compound 
expansions. This is a further development of our scheme for analysis of stability and exponential convergence 
of the PML method [7], see also [9]. The added difficulties are due to the new effect mentioned above. 
To overcome these difficulties we develop our approach to construction of PMLs and use non-homogeneous 
boundary value problems in localization of the essential spectrum and in compound expansions. 

This paper is organized as follows. In Section[5]we introduce manifolds with axial analytic quasicylindrical 
ends and consider an illustrative example. Section[3]is devoted to construction of infinite PMLs. In Section|4] 
we localize the essential spectrum of the operator modeling infinite PMLs conjugated with exponent. In 
Section [5] we establish a limiting absorption principle and show that outgoing and incoming solutions are of 
some exponential decay in PMLs. Finally, in Section [5] we study the problem with finite PMLs and prove 
stability and exponential convergence of the PML method. 



2. Manifolds with axial analytic quasicylindrical ends 

Let S7 be a compact (not necessarily simply connected) n-dimensional manifold with smooth boundary 
dfl. Denote by II the semi-cylinder R + x fl, where M + is the positive semi-axis, and x stands for the 
Cartesian product. Consider an oriented connected n + 1-dimensional manifold M. representable in the form 
M. = M. c Uil, where M. c is a compact manifold with smooth boundary, see Fig.Q] We also assume that the 
boundary dM of M is smooth. 

Let g be a Riemannian metric on M.. We identify the cotangent bundle T*LI with the tensor product 
T*R + eg) T*f2 via the natural isomorphism induced by the product structure on II. Then 

(1) g |"n= Bodx ® dx + 2fli ® dx + g 2 , Bk(x) € C' 00 T*r2® fc , x e E+. 

Denote by CT*fi® fc the tensor power of the complexified cotangent bundle CT*f2 with fibers CT*fi = 
T*fl ® C. In what follows C m stands for sections of complexified bundles. We equip the space C 1 T*f2® fe 
with the norm 

(2) ||-||e=max(|-| e (y) + |V-| t (y)), 

yefi 
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where e is a Riemannian metric on fi, | • \ t (y) is the norm in CT;fi® fc , and V : C 1 T*n Slk -> C*°T*f7® fe+1 is 
the Levi-Civita connection on (Q, e). The norms induced by different metrics e are equivalent. 

Definition 1. We say that (M,g) is a manifold with axial analytic quasicylindrical end (II, g \u), if the 
following conditions hold: 

i. The coefficients R + 3 x i— > Qk(x) £ (7°°x*f2® fe in (fT]) extend by analyticity from the semi-axis R+ to 
the sector S a = {z £ C : | argz| < a < 7r/4}. 

ii. TVie values \\qq(z) — l|| e , ||0i(z)|j e , a71 ^ ||fl2 (^) — f)l!c converge to zero uniformly in z € S a as |z| — » oo, 
where I) is a metric on fl. 

Long-range geometric perturbations of tubular waveguides (i.e. of manifolds (.M, g) with g |~n= da;<g>d:r-|-f)) 
are included into consideration as there are no assumptions on the rate of convergence in condition ii of 
Definition [TJ Let us give some illustrative examples of manifolds with axial analytic quasicylindrical ends. 

Let 11 be a bounded domain in M™ with smooth boundary (£1 is a line segment if n = 1). By (x,y) 
and (s,t) we denote the Cartesian coordinates in R n+1 , where x, s £ R and y,t £ R™. Consider a closed 
domain M. with smooth boundary such that {(x, y) £ M : x ^ 0} is a bounded subset of the half-space 
{(x, y) £ R n+1 : x > -2} and {(x, y) £ M : x > 0} is the semi-cylinder II. Let 

G = {(s, t) £ R" +1 : (s, t) = cf>(x, y), (x, y) £ M} 

be the image of M. under a diffeomorphism <fi satisfying the following conditions: 

i. The function x i-> <f>(x, •) £ C°°(f2) has an analytic continuation from R + to the sector S a ; 

ii. The element in the row i and column m of the matrix (<fi'(z, -))*0'(z, •), where (</>')* is the transpose 
of the Jacobian <fi' , uniformly tends to the Kronecker delta 8i m in the space C 1 (J7) as \z\ tends to 
infinity, z £ § a . 

For instance, we can take 

tf>(x, y) = (x, (x + a + (1 + (x + 3) 7 )?/) , a e M", /3 < 1, 7 < 0. 

Then the boundary dQ approaches at infinity the bent semi-c ylinder {(s,t) : s £ R + ,t - (s + 3) /3 a £ Q}. If 
we take 

4>(x, y) = Q (1 + 1/ log(5 + 4)) da, (l + 1/ log(.x + 5))y^) , 
then 9C/ slowly approaches at infinity the semi-cylinder R + x dft, cf. Fig. [5] Time-harmonic acoustic 




FIGURE 2. Domain £ C M" +1 with quasicylindrical end. 

scattering at angular frequency uj in the domain Q filled by an inhomogeneous medium with anisotropic 
density tensor p and bulk modulus K is modeled by the Neumann Laplacian on (£/, Kp) with spectral 
parameter p = cj 2 , where Kp is a metric on Q, e.g., [3]. Under certain assumptions on regularity of p and 
K the pullback g of Kp by the diffeomorphism <j> satisfies the conditions of Definition [T] These assumptions 
are apriori met for homogeneous isotropic media, i.e. the pullback g of the Euclidean metric Kp by the 
diffeomorphism <f> satisfies the conditions of Definition [T] To study the Neumann Laplacian on (Q, Kp) is 
the same as to study the Neumann Laplacian on g). For our aims it is natural and convenient to work 
on manifolds with axial analytic quasicylindrical ends. 
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3. Construction of infinite PMLs 

Our approach here is very close in spirit to [T3]. Let s r (x) = s(x — r), where r > is sufficiently large and 
s G C°°(K) possesses the properties: 

s(a;) = for all x ^ 1, 

(3) 

^ s'(x) ^ 1 for all and s'(x) = 1 for large x > 0; 

here s' = ds/dx. Let T'*S Q be the holomorphic cotangent bundle {(z,cdz) : z G § a ,c G C}, where 
dz = dHz + id^sz. Consider the tensor field 

(4) g dz ®dz + 2gi <g> dz + 2 2 G C°°(T'*§ Q ® T*n)® 2 

with analytic coefficients S a 9 z i-> Qk(z) G C°°T*fi® , c f Definition [TJ For all values of A in the disk 

(5) O a = {A G C : |A| < sina < l/y/2} 
the function s r defines the embedding 

(6) T*R+ 3 (x,adx) h> (x + Xs r (x),a(l + Xs'^x))' 1 dz) G T'*§ Q , 
where |1 + AsJ.(x)| > 1 — l/\/2. This embedding together with (jJJ induces the tensor field 

(7) g A tn= Qo,xdx ®dx + 2qI x ® dx + g r 2A € C°°T*II® 2 , 

where Q r k x (x) G £'°ot/*£7®' c are sm00 th in x G K.+ and analytic in A G O a coefficients given by 

Qk,x(x) = (1 + Xs' r {x)) 2 - k g k (x + Xs r (x)). 

Note that gA fn with A G O a n 1 is the pullback of the metric g \u by the diffeomorphism x\(x,y) — 
(x + Xsr(x),y) scaling the semi-cylinder LT along its axis R+. Since s r is supported on (r, 00), the equality 
gA t(o,r)xfi= g t(o,r)xn n °l ds for al l A € O a . We extend g A fn to M by setting g A \m\u= g \m\u- As a result 
we obtain the analytic function 

O a 3 X^g x E C°°T*X® 2 . 
Clearly, go = g and (J7J with A = is the same as (fTJ). By analyticity in A we conclude that gA is a 
symmetric tensor field. The Schwarz reflection principle gives gX = gj- It must be stressed that gA with 
A 7^ depends on r, however we do not indicate this for brevity of notations. As r is sufficiently large, gA is 
non-degenerate; i.e. for p G A4 and any nonzero a G CT P M there exists b G CT p A4 such that g^[a, b] 7^ 0, 
where g^[-, •] is the sesquilinear form naturally defined by gA- Indeed, on A4 \ (r, 00) x fl the tensor field gA is 
non-degenerate because it coincides there with the metric g, while on (r, 00) x £1 it is only little different from 
the non-degenerate tensor field (1 + Xs' r ) 2 dx ® dx+ f), see ([7]) and condition ii of Definition [T] In particular, 
gA \n coincides with (1 + Xs' r ) 2 dx (g> dx + f) in the case of manifold with cylindrical end, i.e. in the case 
g |"n= dx® dx + t). 

It is well known that a metric induces the musical isomorphism between tangent and cotangent bundles. 
Similarly, gA induces the fiber isomorphism 

ct;^9^ a kgct p m, 

where A |j£ is a unique vector satisfying the equality £a = g^[ A jj£, a] for all a G CT p M. We extend g^[-, •] to 
the pairs (£, 77) G CT*M x CT;.M by setting 

If A is real, then g^[-, •] is the positive Hermitian form induced by the metric gA. The corresponding volume 
form dvoU extends by analyticity to all A G O a . Introduce the deformed global inner product 

(&w)x= / i\[ZM dvol A , AgO q , Z,ueC?T*M®\ k = 0, 1. 
JM 

Let us stress that for non-real A G O a the form dvolA is complex-valued and the deformed inner product 
(£,cj)a = (o;,^)a is not Hermitian. Let us consider the sesquilinear quadratic form 

qx[u,u] = (du,d(g % u)) X7 ugC c °°(X), A G O a , 

where d : C^°(M) — > C^°T*M is the exterior derivative and g\ G C°°(M) is such that g\ dvoU = dvolo. 
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Lemma 1. Introduce the space L 2 (Ai) and the Sobolev space H 1 (A4) as completion of the set C^°(A4) with 
respect to the norms 

\\u\\ = y/(u, u)o, \\u\\ H i( M) = ^J(du,du)o + (u,u) 

respectively. Then the family O a 3 A i— > q\ of unbounded quadratic forms in L 2 (A4) with domain H l (Jv[) 
is analytic in the sense of Kato [101 116) : i.e. is a closed densely defined sectorial form, and the function 
O a 3 )h> c\x[u, u] is analytic for any u € H 1 (A4). Moreover, the sector of c\\ is independent of A £ O a . 

The family O a 3 A H qj uniquely determines an analytic family O a 3 A n- A A of unbounded m- 
sectorial operators in L 2 (M.) [10] [16]. (Here and elsewhere m-sectorial means that the numerical range 
{(Au, u) : u € @(A)} and the spectrum a(A) of a closed unbounded operator A with the domain &(A) are 
both in some sector {/x € C : | arg(/i + c)| ^ < tt/2, c > 0}.) We have 

(8) ( X Au,v) =qx[u,v], ue@{ x A),veH\M). 

In particular, (°Au,v)o = (du,dv)o and °A is the selfadjoint Neumann Laplacian. We consider A A with 
non-real A € O a as the operator modeling an infinite PML on (r, oo) x f2. We will show that this PML is an 
artificial nonreflective strongly absorbing layer for the outgoing (resp. incoming) solutions if 3A > (resp. 
QA < 0). Note that for u G f^(°A) supported outside (r, oo) x we have u G 3>( X A) and °Au = x Au. 

Proof. In a finite covering of fi by coordinate neighborhoods {^j} take a neighborhood ^ and coordinates 
y e 1" in Then {IR + x ^ } is a covering of II and (x,y) are coordinates in M + x ^ . Denote do = 
d/dx, di — d/dye, do — dx, and de = dyg. As gA is non-degenerate, the matrix g\/ m in the coordinate 
representation gA = g\,& m di ® d m has the inverse g A m and for £ = we have A j}£ = g£ m &9m- Hence 

gAi'.;- gri-c,,. 

On f2 all metrics are equivalent. In particular we can take e = di m de ® ci m , where <5f m is the Kronecker 
delta. Then Definition Q] together with (J7J immediately implies that 

(9) |^(gi m -diag{(l + A)- 2 ,r 1 }£ m )(x, 2 /)| C(x) ^ as x ^ oo, j = 0, 1. 
Moreover, 

(10) |^(gi m -diag{(l + A4)- 2 ,r 1 W)(x,j/)| <e,, x^r, A e O a , j = 0, 1, 

with sufficiently small e r as r is sufficiently large. Since f) (y) is a symmetric positive definite matrix, we 
get 

|arg((l + A S ' r (x))- 2 |Co| 2 + t) lm (y)te m ) \<2a< tt/2, 

ciei 2 < |(1 + A4(x))- 2 |e | 2 + f) /m (j/)6£ m | < \e/c, 
where |£| 2 = £^ and c > 0. This together with flTU|) gives 

(11) ^ICI 2 < %A m (*,?/fe) < <^| 2 , I arg(g £ A m (x,y)^" m )l < 0, 

where x ^ r, S = min{c — e r , (1/c + e,,)^ 1 }, and d = 2a + 2 arcsin(e r /2c) < 7r/2. 

From (ITT1) and g A = g p , p € .M \ (r, oo) x il, we conclude that the first term in the representation 

(12) q A [tt,u]= / g A [du,du]dvolo+ / i^l^ldvolo, ueC^(M), 

J m J m Pa 

meets the estimates 

(13) — 9 ^ arg / g\[du, du] dvolo ^ 0, c{du,du)o ^dt g\ [du, du] dvolo (du,du)o/c 

Jm Jm 

with some 9 < n/2 and c > 0. The second term has an arbitrary small relative bound with respect to the 
first term. Indeed, 



(14) 



I g x [du,ud^ dvolQ ^ (1+tang) / /• gA[ ^ ]dvol V /2 

X ($t [ \u\ 2 g\[dQ\,dQ\] dvolo ^ s% e / gA[du, du] dvol 
V / ./ai 



+ £- 1 C||m|| 



for any e > and an independent of e, A, and u constant C. Here we used the uniform in A and p estimates 
|f?A(p)| ^ c > and g^[dg\,dg\] ^ C < oo, which are valid because Q\(p) — 1 for p £ (r, oo) x f2 and 
Qx = \/\go\/\g\\ in R+ x % % where d\\fc > \{x,y) -> + Xs' r ) 2 H)\)(x,y) as x > r -> oo, cf. (JTO]) (as usual, 

|g| := det(g fm )). 

As a consequence of (| 1 2[) . (|13[) . and fj 14|) for some i) < tt/2 and some positive S and 7 we obtain 

(15) |arg(q A [u,M] +7||w|| 2 )| < d, S(du, du) - -f\\u\\ 2 ^ ftqx[u,u] ^ ((du, du) Q + \\u\\ 2 )/S 

uniformly in A £ O a . Therefore q\ in L?(M) with domain H 1 (Ai) is a closed densely defined sectorial form 
and its sector {fi £ C : | axg(fi+j)\ ^ #} is independent of A. By construction the function O a 9 A >->• q\[u, u] 
is analytic for any u £ H 1 (A4). □ 

4. Conjugated operator and its essential spectrum 

In order to show exponential decay of outgoing (resp. incoming) solutions in infinite PMLs we study the 
operator A A, 3A > (resp. 3A < 0), conjugated with an exponent. 

Let s be a smooth function on the semi-cylinder II, which depends only on the axial variable x £ R+ and 
possesses the properties (j3|). We extend s to a smooth function on M. by setting s tM\n= 0. Consider the 
conjugated operator A A^ = e^ sA Ae fls with parameter (3 £ C on functions {u £ C^°(M) : e@ s u £ ^( A A)}, 
where e@ s is the operator of multiplication by the exponent. With x Ap we associate the quadratic form 



qf [u,«] = (d(e^ s u),d(e-^ A u)) A) u £ C™(M). 

Lemma 2. For any A £ O a the form q^ in L 2 (M) is a closed sectorial form with the domain H 1 (M.). 
Moreover, its sector is independent of X. 

Proof. We have 

q A [u, u] — c\\[u, u] = f3(uds, d(gju)) . — (3(g\ du, uds) — (3 2 (g\uds, uds) x , 

where the right hand side depends linearly on du. Similarly to (fl"4"]l we conclude that the difference — q\ 
has an arbitrarily small relative bound with respect to q\. More precisely, 

|qf [u,u] - qx[u,u}\ < e|qA[t*,u]| + C(\f3\, e)\\u\\ 2 , 

where e > is arbitrarily small and C(|/3|,e) is independent of A £ O a and u £ C ( ?°(A / J). This together 
with (fT5"]) completes the proof. □ 

The Friedrichs extension of A A^ is an m-sectorial operator [10] [16]. Consider its domain S>( x Ap) as a 
Hilbert space with the norm • || 2 + || A A. j a • || 2 . We say that p is a point of the essential spectrum (r es8 ( x Ap) 
if the bounded operator 

(16) x Ap - M : @( X Ap) -> L 2 (M) 

is not Fredholm. (Recall that a bounded linear operator is said to be Fredholm, if its kernel and cokernel 
are finite-dimensional, and the range is closed.) 

Proposition 1. Let A £ O a and (3 £ C. Then u £ a- ess ( Ap) if and only if 

(17) Ai -(l + A)- 2 (e + z/3) 2 e ( x(A^) 

for some ^el, where <t(Aq ) is the spectrum of the selfadjoint Neumann Laplacian Aq on (ft, f)). 

The spectrum cr ess ( A A^) is depicted on Fig. [3J In the case (3 = the parabolas collapse to the dashed 
rays and we obtain the essential spectrum of A A = A Ao . The proof of Proposition [T] is preceded by 

Lemma 3. The operator A A, A £ O a , corresponds to a regular elliptic boundary value problem on A4. 

Proof. On A4 \ (r, 00) x SI the operator A A coincides with the Neumann Laplacian, which corresponds to a 
regular elliptic problem [15] [13] . 
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FIGURE 3. Essential spectrum of the operator A A,g for 3 A > and (3 ^ 0, the dashed rays 
correspond to /3 = 0. 

Let y G R" be a system of coordinates in a neighborhood ^ on fi. In the case dft fl f j / we pick 
boundary normal coordinates y = (y',y n ) such that d Vn coincides with the unit inward normal derivative 
given by the metric t). From ([5]) it follows that in the coordinates (x,y) we have 

(is) a a - -igAr i/2 ^| gA r/ 2 g i m c? m 

and the functions u € £F( A A) satisfy some boundary condition x Nu = on dAi such that 

(19) x Mu = g n x m d m u \ Vn = = 

if dil C\ a i/j 7^ 0; here notations are the same as in the proof of Lemma [JJ The operator (TIB)) is strongly elliptic 
for x ^ r due to It remains to show that the Shapiro-Lopatinskh condition is met on (r, oo) x (dflD^j). 
Thanks to (fTQ|) the principal part of ([TBI) f° r Vn = (resp. the principal part of the operator of boundary 
conditions in (Tl9|) ) is little different from —(1 + \s' r (x))~ 2 dQ — Q(y',d y ') — d 2 (resp. from d n ) uniformly in 
x ^ r, A G 2? Q , and y' as r is large. Here —Q(y\d y ') — d% is the principal part of An- For £ = (£o,£') on 
the unit sphere S n and every bounded solution u ^ of 

((1 + As;(x))- 2 £ 2 + Qfj/,£') - d 2 Ju(y n ) = for y„ G K+ 

we have 

d Vn u(0) = Cd yn e-y^(i+^))- 2 e B +Q(y'^) ^ =0= + As^ar))-^ 2 + Q(j/', f) ^ 

because 5R(1 + As^.(cc))~ 2 > and Q(y' ,£,') 1S a positive dehnite quadratic form in In other words, the 
pair (—(1 + \s' r (x))~ 2 d 2 + Q(y',d y >) — <9 2 ,<9„) satisfies the Shapiro-Lopatinskh condition or, equivalently, 
the estimate 

(20) l|u|| ff2(R+) < C(||((l + \s' r (x))- 2 e a + Q(y',e) - dl)u\\ L2{R+) + \d yn u(0)\) 

holds, e.g. [131 US]- The constant C in ([20]) is independent of r > 0, £ G 5™, A G O^, x G K + , and 
y 1 G ^ H Sf2. (Indeed, it suffices to note that (|20p with r replaced by r' can be obtained by the change of 
variables x h-> x+r— r' , the function s£ varies over a compact subset of R+ only, and (l+Xs' r (x))~ 2 £ l Q+Q(y' , £') 
is smooth in A, £, a;, and y'.) This together with ([TU]) implies that the estimate similar to ([2TJ)l is valid for the 
principal parts of dTBl and (TK))) provided x ^ r, where r is sufficiently large. This completes the proof. □ 

Now we are in position to prove Proposition [T] 

Proof. We will rely on the following lemma due to Peetre, see e.g. [T51 Lemma 5.1] or [T31 Lemma 3.4.1]: 

Let X,y and Z be Banach spaces, where X is compactly embedded into Z. Furthermore, let £ be a 
linear continuous operator from X to y. Then the next two assertions are equivalent: (i) the range 
of C is closed in y and dimker£ < oo 7 (ii) there exists a constant C such that 

(21) \\u\\ x <:C{\\Cu\\ y + \\u\\ z ) \fueX. 

Sufficiency. Here we assume that /i does not satisfy flTTl) and establish an estimate of type ([2"Tj) for the 
operator OH). Consider the operator ( A A, A W) : C™(M) -> C%°(M) X C™(dM) satisfying 

{ x Au,v) + ( x Afu,v) =qx[u,v], ueC^°(M),v£H 1 (M), 
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where W is the operator of boundary conditions on dM and (•, •) is the inner product on {dM.,g \dM.)- 
Clearly, x Nu — for u £ ^( A A). Besides, in the local coordinates (x,y) on M + x we have (fT8]) and (fT9j) . 
In particular, °A is the Laplace-Beltrami operator on (.M,g) and IV is the corresponding operator of the 
Neumann boundary conditions. Consider also the conjugated operator ( A A^, "Wa) = e _5/3 ( A A, 'W)e s/3 . 

Introduce the Sobolev space H £ (M. x fi) as the completion of the set C^°(M. x fi) with respect to the norm 

IMIff'(Rxn) = ( / ^II^M^II^-'-cn)^) 7 > 

where H l (Vl) is the Sobolev space on n. Let jZV2(R x <9f2) be the space of traces u tRxdn of the functions 
u G i? 1 (R x f2). Denote by the operator of Neumann boundary conditions on M x dVt taken with respect 
to the metric dx ® dx + f) . Applying the Fourier transform J- x ^ we pass from the continuous operator 

(22) (A n - (1 + X)~ 2 (d x + /3) 2 - n,d n ) : H 2 (R x ft) -> ff°(R x ft) x H 1/2 (R x 90) 

to the operator (An + (1 + A)~ 2 (£ + i/3) 2 — (J,, d v ) of the Neumann problem on (ft, f)). As /j, does not meet (|17p . 
the latter operator is invertible for all £ and the inverse of (|22[) is given by 



^(An + (1 + A)- 2 (e + iP? - /i, 9,) Y T, 



see e.g. [131 Theorem 5.2.2] or [T2J Theorem 2.4.1] for details. As a consequence we have 

(23) IMU'CRxn) < C|| (An - (1 + A)" 2 ^ + p) 2 - fx, d v )u\\ H0{Mxn)xHl/2(Kxgay 

Let Xt(x) = x( x ~~ T), where \ <= C°°(R) is a cutoff function such that x( x ) = 1 f° r x ^ 1 and x(x) = 
for x < 0. Then due to ©, iflg]). and ([19]) the constant c(T) in the estimate 

|| ( A A^ - An + (1 + \y 2 (d x + I3) 2 , x Afp - a i? )xTu|| H o (Rx n)xHi/2 (K xao) < c(T)\\xtu\\ H ^(Rxn) 
tends to zero as T ~ > +oo. This together with (|2"31) implies that for all sufficiently large T the estimate 

(24) \\XTu\\m(MxQ) ^ C \\( X A/3 - Ll, X J\ffj)XTu\\H«(Rxfl)xm/2(Rxdn) 

holds, where C = (1/C - c(T))" 1 > 0. 

Let 7J f (A4) be the Sobolev space introduced as the completion of the set C%°(M) in the norm 

\\ u \\h*(m) = (\\u\\ 2 H t {Mc} + I ^\\^{x)\\ 2 H i- k(sl) dx^ ' , 

where d x — d/dx and H e (Ai c ) is the Sobolev space on the compact manifold A4 C - By i/ 1 / 2 (<9A / f) we denote 
the space of traces on dA4 of functions in H 1 (A4). We extend xt from its support in LI to M. by zero and 
rewrite (l24l) in the form 



\\xtu\\ H 2 (M ) < c(||xt( A A,3 - ^, A A/>)u|| L 2 (A4)xH i/2 (A1) + ||[( A A,g, x Afp),XT}u\\ L 2 {M)xH i/2 (M) . 
For the commutator we have [( A^, x Np), xt]u = [( A A,3, A A/]a), Xt](1 — X2t)m and therefore 

|| [( A A^, A A/»,Xt]m||l2(X)xH1/2(> < ) < C||(l - X2t)u\\ H 2(M)- 

Moreover, as a consequence of Lemma [3] the local elliptic coercive estimate 

||(1 - X2t)u\\ H 2(M) < C(l (1 - X3t)( A A,3 - (I, A A/»M|| L 2 (A4)xff i/2 (aA4) + 11(1 - X3t)w||) 

is valid [15]. From the last three estimates it follows that 

(25) \\u\\ H 2(M) < C(ll( A A/9 - H, X Ni3)u\\ L 2 iM)xH i/2 {dM) + 11(1 - X3t)w||)- 

In particular, ([25]) implies that || • ||_ff 2 (x) is an equivalent norm in ^( A A^) = {ii£ H 2 (A4) : x Mpu = 0}. 

Let w be a bounded rapidly decreasing at infinity positive function on Ai such that the embedding of 
H 2 (M) into the weighted space L 2 (A4,w) with the norm ||w • || is compact. Then ^( A A^) is compactly 
embedded into L 2 (A4,w) and ([25]) is an estimate of type (|2"Tj) for the operator ([i~6]) . Thus the range of (fi"6]) 
is closed and the kernel is finite-dimensional. In order to see that the cokernel is finite-dimensional, one can 
apply the same argument and obtain an estimate of type (|2ip for the adjoint operator. 

Necessity. Now we assume that \x meets (|17[) for some j and show that the operator (|16|) is not Fredholm. 
It suffices to find a sequence {ve}f^ 1 in S>( x /S) violating the estimate ([25]) . 



We first show that for a regular point /j,o of the m-sectorial operator A A^ the continuous operator 

(26) ( A A,g - u , A A/» : H 2 (M) -> L 2 (M) x H 1/2 (dM) 

realizes an isomorphism. With this aim in mind we replace (1 — xst)u by wu in (|25|) . Then by the Peetre 
lemma the range of (|26p is closed. It is easy to see that the elements in the cokernel of the operator (|2^|) 
are of the form (v,v Iom) with v € ker( A A^ — (iq), where A A^ is adjoint to the m-sectorial operator x Ap in 
L 2 (A4) with domain £^( A A^). Indeed, let (v,v_) be in the kernel of the adjoint operator 

( A A,g — fiQ, x J\fp)* : L 2 {M) x {HV\dM)Y (H 2 (M))*. 

Then ( x Apu - fi u,v) + ( X Npu,v) = for all u € H 2 {M), where (•,•) is extended to H 1/2 (dM) x 
(ifV^At))*. Since ^( A A/3) C H 2 (M), we immediately see that u G ker( A A^ - u ). Then for u £ H 2 (M) 
the Green identity gives ( x Afi3U,v — v \om) = and therefore v — v \om- As a consequence, (|26|) is an 
isomorphism for /xo ^ <r( A A^). 

Let x € C°°(R) be such that x(ir) = 1 for |i — 3| < 1, and x( x ) = for |a; — 3| ^2. Consider the functions 

(27) ui(x, y) = x(x/£) exp(i(l + \)y/fi - Vjx)&(y), (x,y)eRxQ, 

where $ is an eigenfunction of the Neumann Laplacian Aq corresponding to the eigenvalue vj . It is clear 
that ui satisfies the Neumann boundary condition d v ug — on R x d£l. As u meets the condition (fTT]) . the 
exponent in ([27]) is an oscillating function of x G R. Straightforward calculation shows that 



(28) ||(A f2 -(l + A)- 2 (^+/3) 2 — M) u ^||jjo(n x n) ^ const, \\ue\\H 2 (Ux(i) 00 
as £ — > +oo. We extend the functions ut from II to Ai by zero and set 

vt =u t - { X A P - Mo , ^^(O, x JV/jU/), 
where is a regular point of the m-sectorial operator x Ap. Clearly, vg £ Si^Ap). We also have 

(29) H^Ap - x ^)- 1 (a, : V^«^)||ji a (A4) < C^IK^ - ^)«elUv a( ajw) < C,|M|„ 2(Rxn) , 
where Cg — > as I — > +oo. Hence 

(30) \\ve\\H 2 (M) ^ ||^||/f 2 (Rxn) - CeWuiWnnRxn), \\vt - ue\\ H 2(M) ^ CU||w*||jy2(R X n)- 

Assume that the estimate ([23)) is valid. Without loss of generality we can take a rapidly decreasing weight 
w such that ||wu^| ^ Const uniformly in I and the embedding H 2 (M.) L 2 (A4;w) is compact. It is clear 
that 1 1 wu 1 1 ^ c||u|| ff 2(_ A/1 -) with some independent of u £ H 2 (M.) constant c. Therefore 

(31) l|wu*|| < Const + cC/||^||//2(Rxn), 
cf. dSTI). The estimate 

||( A A,3 - A n + (1 + \y 2 (d x + pf)u t \\ ^ c/Mfla^n), 
where eg — > as i — > +oo, together with ([2"5)l and ([317)) gives 

|| ( A A^ - M M| < C\\ (-(1 + X)- 2 (d x + p) 2 + An - mHIho (R xO) 

(32) + ||( A A^ - An + (1 + \)- 2 (d x + p) 2 )u f \\ + \\{ X A - n){ Vi - u e )\\ 

^ C ■ const + (a + CC e )\\ui\\ H 2( Rxfl ). 

Finally, as a consequence of ([30]). ([25]) . ([32]) and (f3Tj) . we get 

^ 33 ^ II^IIh 2 (rxI2) - C^||^|| H 3(Kxn) < \\vi\\ H 2(M) < C(||( A A^ - || + ||wv £ ||) 

C(C • cortsi + (CQ + c/)||ti/|| H 2 (Rxn) + Const + cC e \\ui\\ H 2( Rxn) ). 



Since q — > and — > 0, the inequalities (j33[) imply that the value (Kxfi) remains bounded as 

I — > +oo. This contradicts ([28]) . Thus the sequence {vi} ( ^ 1 violates the estimate ([25]) . □ 
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5. Exponential decay of solutions in infinite PMLs 



Consider the algebra & of all entire functions C3zh> F(z, •) £ C°°(f2) with the following property: in 
any sector \$sz\ ^ (1 — e)3?z with e > the value \\F(z, -)llL 2 (f2) decays faster than any inverse power of Sftz as 
dtz — » +oo. Examples of functions F € § are F(z, y) = e~ Jz2 P(z, y), where 7 > and P(z, y) is an arbitrary 
polynomial in z with coefficients in C°°(f2). We say that / £ L 2 (A4) is an analytic vector, if f(x, y) = F(x, y) 
for some F £ S and all (x, y) £ IT. For any / in the set A of all analytic vectors we can define the analytic 
function O a 3 X h-> fx £ L 2 (M) by setting fx = f on M \ (r, 00) x fi and fx(x,y) = f(x + Xs r (x),y) for 
(x,y) £ (r, 00) x fi. The set {fx : / € -4} is dense in L 2 {M) for any A £ T> a , see, e.g., [51 Theorem 3]. 

Theorem 1. Assume that /to £ R\er(Aj}) is not an eigenvalue (i.e. not a trapped mode) of the self-adjoint 
Neumann Laplacian °A in L 2 (A4). Then the following assertions are valid. 

1. For any f £ A there exist outgoing u_ and incoming u + solutions defined by the limiting absorption 
principle 

(34) u+ = lim(°A - fi - ie)- x f, u_ - lim(°A - Mo - ie^f, 

where the limits are taken in the space Lf oc (A4). 
2. For A £ O a with 3A < (resp. 3A > 0) the m-sectorial operator A A in L 2 (A4) models PMLs on 
(r, 00) x f2 /or incoming (resp. outgoing) solutions in the sense that ux — ( A A — l^o) 1 f\ coincides 
on M\ (r, 00) x 17 wii/i t/ie incoming solution u + (resp. with the outgoing solution U-) and 

(35) \\e 0s ux\\ H ^ M ) <C(^,X)\\e^fx\\, < /3 < rmn uea(AS) \^{(l + A)V^U}|. 

TTie estimate (|35[) shows that ux decays exponentially in the PMLs. 

Proof. 1. We need to prove that for any £ € C%°(M) and / e A the function g(°A — (iq — ie) -1 / tends to 
some limits in L 2 (A4) as e 4- and e f 0. Take a sufficiently large r such that supp g n (r - , 00) x f2 = 0. 
Let /x be outside the sector of m-sectorial operator A A. Then for any real A £ O a the change of variable 
x 1 y x ~\~ As r (a;) outside supp g implies 

(36) e(°A- M )- 1 / = p( A A-/i)- 1 /A. 

From Lemma[T]and ([8]) it follows that the resolvent ( A A — /i) _1 is an analytic function of A £ O a with values 
in the space of bounded operators in L 2 (Ai) [T6J Theorem XII. 7]. Thus (|36|) extends by analyticity to all 
A £ O a - As /i is a regular point of A A in the the simply connected set C \ a ess ( x A) (see Proposition [T]) , the 
Fredholm analytic theory implies that the resolvent 

C \ a ess { x A) ( A A — ^ : L 2 (M) -> ^( A A) 

is a meromorphic operator function; see, e.g., |12[ Appendix]. It remains to show that no £ C \ cr ess ( x A) is 
not a pole. Then the right hand side of ([55]) with 3A < (resp. 3A > 0) provides the left hand side with 
analytic continuation from /x = /to + ie to /i = /xq as e I (resp. e f 0). 

For /x outside of the sector of A A and real A £ O a by the change of variable x H> x + Xs r (x) we obtain 

(37) {(°A-n)-if,g) = (fA-n)- 1 fx,9x)x> f^eA. 

This equality extends by analyticity to all A £ O a . Suppose, by contradiction, that /io is a pole of ( A A — /i) _1 
with 3 A sg 0. As the sets {fx-fEA} and {gx ■ g £ A} are dense in L 2 (M), the right hand side of ([37]) 
has a pole at /xo for some / and g. Then (|37|) implies that (P^ /, 5) ^ 0, where P Mo is the projection onto 
the eigenspace of the selfadjoint operator °A, and thus ker(A — /xo) 7^ {0}. This is a contradiction. 

2. Since (|3"rJj) is valid for any g £ C%°(M.) such that supp g n (r, 00) x f2 = 0, our construction in 
the proof of assertion 1 shows that ux coincides on M. \ (r,oo) x with u+ (resp. u_) if 3A < (resp. 
3 A > 0). The condition on (3 in (|35|) guarantees that /io together with all points outside of the sector of 
m-sectorial operator x Ap is in the simply connected subset of C \ a ess ( x Ap), see Proposition [TJ Then the 
Fredholm analytic theory implies that /to is either an eigenvalue oi x Ap or /i ^ a( x Ap), e.g., [12j Appendix]. 
The inclusion "J/ £ ker( A A ( g — /x ) gives e _/3s \I/ € ker( A A — /t ) = {0} as /t is not a pole of ( A A — /u) -1 . 
Hence /to ^ cr( A A / 3). This together with the equality e" s itA = { x Ap — no)~ 1 e fx justifies the estimate ()35|) . 
cf. PSJ|. □ 

in 



6. Finite PMLs, stability and exponential convergence of the PML method 

Consider the compact manifold Mr = M \ (R, oo) x 57. The boundary dM r oIMr is piecewise smooth. 
It has two conic points in the case of a 1-dimensional manifold 57 and the edge 957 x {R} otherwise. We 
denote OMr \ ({R} xfl) = dM^- Introduce the Sobolev space H 2 (Mr) as the completion of the set 
C£°(Mr) with respect to the norm 

\\u\\h2(M r ) = (\M\h2(m c ) + 12 / W d i u ( x )W 2 Hi-*(n) dx ) ' 

where H 2 (M C ) and H 2 ~ e (Q) are the Sobolev spaces on the smooth compact manifolds M c and 57. Consider 
the problem with finite PMLs: given F £ L 2 (Mr) find a solution v G H 2 (Mr) to the problem 

(38) ( A A — no)v = Fon -Mr; x Nv = on dM R ; v = on {R} x 57. 
In the next theorem we prove a stability result for this problem. 

Theorem 2. Let X € O a \ R. Assume that /iq £ M \ cr(An) is not an eigenvalue of the selfadjoint Neumann 
Laplacian °A in L 2 (M). Then there exists a large number Rq > such that for all R > Rg and F G 
L 2 (Mr), (|38p has a unique solution v £ H 2 (Mr). Moreover, the estimate 

(39) \\v\\ H 2 (Mr) ^C\\F\\l^m r ) 
holds with independent of R > Rq and F constant C . 

In the proof of Theorem[2]we rely on compound expansions. This requires construction of an approximate 
solution to (|3"51) compounded of solutions to limit problems, e.g., [T3]- Being substituted into (|3"5|) the 
approximate solution leaves a discrepancy, which tends to zero as R increases. In contrast to the case of 
the Dirichlet Laplacian [7], here the discrepancy left in the boundary conditions on dM^ cannot be made 
small for large R if we use homogeneous limit problems. As the first limit problem we take the problem with 
infinite PMLs and non- homogeneous boundary conditions on dM. In the next lemma we study the second 
limit problem. 

Introduce the weighted Sobolev space oo, R) x 57) as the completion of the set C^°((— oo, R) x 57) 

in the norm 

1/2 




H*((-oo,ii)xn) = 2^/ H e d x^x)\\ Hk - e{n) dx 



We also set L 2 ((— oo,i?) X 57) = H^((—oo,R) x 57) and denote by H^ 2 ((—oo,R) x 957) the space of traces 
on (— oo,R) x 957 of the functions in Hi((— oc,R) x 57). 

Lemma 4. Assume that A <E Oq, \ R, /io £ 1 \ <j(Aq), and /3 is in the interval (|35|). XTien /or any 
f € oo, -R) x 51) and g € H^ 2 ((—oo, R) x 957) £/iere exists a solution u to the second limit problem 

(40) (A n - (1 + A)~ 2 9 2 - Ho)u =f on (-oo,R) x 57, 9,,u = g on (-oo, i?) x 957, u = on {i?} x 57, 
satisfying the estimate 

( 41 ) || u llff|((-oo„R)xfi) < C r (l|f||i|((-oo,ii)xfi) + Hgll//i/ 2 ((-oo,_R)xc3n))' 

where the constant C is independent off and g. Here d v is the operator of the Neumann boundary conditions 
on (— oo, R) x 957 induced by the product metric dx ® dx + f) on (— oo, i?) x 57. 

Proof. Without loss of generality we can assume that R — (the general case can be obtained by the change 
of variables x h-> x — R). We give only a sketch of the proof as it is essentially based on a well known 
argument, details can be found e.g. in [13] . 

Assume that (f , g) e C7?° (R_ x 57) x C™° (R_ x 957) and extend it to a function in C c °° (R x 57) x C c °° (R x 957) 
by setting (f, g)(— x) = — (f, g)(x) for a; < 0. The Fourier transform T xi ^ £ (f , g) is an entire function of £ 
with values in L 2 (57) x iJ x / 2 (957) decaying faster than with any fc > as £ tends at infinity in any 

strip < C. Since /io — (1 + A) _2 £ 2 with ^ 3£ < f3 is not an eigenvalue of Aq and the strongly elliptic 



operator An — (1 + A) 2 <9 2 on R x with the operator d n on R x d£l set up a regular elliptic problem, the 
elliptic coercive estimate 
2 

(42) l£l 2 l*H^-(o) < C(ll((l + A)- 2 C 2 + Ao - n)nl H u) + \t\\\W\lnoa> + \\W\%t/ H oo)) 

with parameter ^ Q£ < (3 is valid. Moreover, the operator 

(A n + (1 + A)- 2 £ 2 - Mo , 9,) : H 2 {Q) -> tf°(0) x H 1 / 2 ^), < 9£ < A 
yields an isomorphism and 

(43) u(x) = 7^ (An + (1 + A)- 2 e 2 - Mo, ^) _1 ^(f, g) 
is a unique in L 2 (M x Q) solution to the Neumann problem 

(Aq - (1 + A)~ 2 <9 2 - ^o)u = f on R x Q, 9,u = gonIx 90. 

Usual argument on smoothness of solutions to elliptic problems gives u £ C°°(M x fi). Since x H> (f, g)(x) 
is odd it follows that T xh ^^(f, g) is an odd function of £. Therefore x i-> u(x) is odd and u(0) = 0. 

The Cauchy integral theorem and the estimate (|42p allow to replace the contour of integration R in 
by {£ G C : 3£ = j3}. Then the Parseval equality implies that u is a solution to the problem (14TJ1) satisfying 
the estimate (|4"Tj) . By continuity our construction extends to all f £ L 2 (R_ x il) and g £ _ff^ 2 (R_ x 951). □ 

Now we are in position to prove Theorem [2] 

Proof. The proof is carried out using the compound expansion method. We say that w £ H 2 (Mr) is an 
approximate solution to the non-homogeneous problem 

(44) ( A A - (i )v — F on Mr, x Mv = G on dM R , v = on {R} x fl 
if the following conditions are satisfied: 

i. The estimate 

II^IU 2 (Mji) < c \\( F i G )\\ L 2(M R )xHi/2(dM R ) 

holds with an independent of F, G, and R constant c; 

ii. The estimate 

(45) ||( A A - ^o^AA)^ - (^,a)|| £2(A4A)xJyl / 2(a ^i-) < C^ll(^^)]li 2 (^)x^/-(a^f-) 
is valid, where the constant Cr is independent of F and G and Cr —¥ as R — > +oo. 

Due to Condition i wr continuously depends on F and G. Condition ii implies that the discrepancy, left by 
wr in the problem (|44p . tends to zero as R — > +oo. Once an approximate solution is found, it is not hard 
to verify the assertion of the theorem. 

Let p £ C°°(R) be a cutoff function such that p(x) — 1 for x ^ and x( x ) = f° r x 1/2- We set 
Pr = p(x — R) on II and qr = 1 on M \ II. Let (/, g) — qr/ 2 (F,G) and (f, g) = (1 - qr)(F,G). Extend 
(F, G) from Mr x dM R toMxdM and (f , g) from (0, R)xflx (0, R) x dQ to (-co, R) X O X (-co, i?) x 90 
by zero. We already know that po ^ cr( A A^) and hence ([2TH) yields an isomorphism for /3 in the interval 
in (1551 (see the proof of Theorem[T]2). We find an approximate solution w compounded of u\ — e _,3s ( A A^ — 
po, x Np)^ 1 e^ s {f,g) and a solution u £ H 2 {{~ oo,R) x Q) to the equation (|4T7|) in the form 

w = qru x + (1 - to/ 3 )u; 

here the second term in the right hand side is extended from (— oo, R) x Vt to Mr by zero. On the support 
of (/, g) we have e^ s ^ Ce^ R / 2 and on the support of (f, g) we have e~^ x ^ Ce~P R l 2 uniformly in R. Hence 

l|e /3S (/,g)||L2(A4)xHi/2 (9M ) + e ' 3i? ||( f 'g)llH«((-oc,fl)xa)x^ /2 ((-oo, J R)xan) 

< Ce fm/2 \\(F, G)\\ L 2 {MR)xH i/2 (dM -- ) 
with an independent of R and (F, G) constant C. Thanks to (14"TT) and (|4"6]) we can conclude that 

\\ W \\h 2 (Mr) ^ II^ U a|Ih 2 (M) + IK 1 ~ PR/3) U \\H 2 ((- o,R)xn) < C 'll( F ' G )llL2 (A4R)xff i/2 (a7WH) 
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and Condition i is satisfied. 

Let us verify Condition ii. We have 

, s ( x A- f i , x Af)w-(F,G) = {( x A, x Af),Q R }u x +((l + \)- 2 {dlp R/3 }u,0) 

(47) 

+ ( A A - An + (1 + \y 2 dl X N - d v ){\ - p R/3 )u. 

On the support of the commutator [( A A, x Af), Qr]u\ the weight e^ s is bounded from below by ce@ R uniformly 
in R > 0. As a consequence, 

(48) || [( X A, x Af), Qr] u \\\ L 2( Mr ) xH i/2( 9M -> ) *S Qe-WWeP'uxW&M < C 2 e'^ R \\e^(f, g)\\ L2(M)xH1/HdM) . 

Now we estimate the second term in the right hand side of (j47|) . On the support of [d 2 ,p R / 3 ]u we have 
e-P x ^ Ce-' 3 ^/ 3 . Thanks to gO we obtain 

\\(l + \)- 2 [dlp R/3 ]u\\ L2(MR) ^C 1 e^ R / 3 \\u;Hj((-oo,R)xn)\\ 

< r BR/3\\/c \|| 

Finally, consider the last term in the right hand side of (14"T1) . On the support of (1 — Pr/ 3 )u the coefficients 
of the operator ( A A - Aq + (1 + A)~ 2 <9^, x Af-d v ) uniformly tend to zero as R ->• +oo; sec © and (fl8|). ([19]) . 
This together with (l4"Tj) gives 

II (" A " A n + (1 + >-r 2 ^, W - 9„)(1 - Pr/M\l W ,h^ { 8M- r ) 

^ c fl||(f'6)llij2((_oo,.R)xn)xffo /2 ((-oo R y Kdn y 

where c R — > as i? — > +oo. From (|4^1) -([5U | it follows that w meets Condition ii. Thus w is an approximate 
solution to the problem (|44| . 

Now we are in position to prove the assertion of the theorem. Observe that ( A A, x Af)w — (/, g) = D R (f, g) 
with some operator D R in L 2 (Ai R ) x H 1 ^ 2 (dA4] i ), whose norm ||D.r|| tends to zero as R — > +oo because of 
Condition ii on w. For all R > Rq we have \\O r \\ ^ e < 1. Hence there exists the inverse (/ 4- Dij) -1 and its 
norm is bounded by the constant 1/(1 — e) uniformly in R > Rq. We set (/,<?) = (I + D R )~ 1 (f,g). In the 
same way as before we construct the approximate solution w for the problem (|4"4"]l . where (/, g) is replaced 
by (f,§)- Then for v = w we have ( A A - p , A A/> = (/,<?) + R {f 7 g) = (f,g) and 

\\v;H 2 (Mr)\\ < c\\(f,9)\\mM R )xm/2(dM R ) < c /( 1 - e )IK/'.9)llL2 (A4fl)xH i/2 (aA4ii) , 

where C is independent oi R > Rq. In particular for R > Rq and / € L 2 (AAr) there exists a solution 
i> G H 2 (M.r) to the problem (1551) satisfying (j3"9")l . This solution is unique as a similar argument shows that 
the adjoint problem is solvable in H 2 (A4r) for any right hand side in L 2 (A4r). □ 

In contrast to infinite PMLs, finite PMLs are not non-reflective. Reflections produce a non-zero difference 
in Ai r between solutions u± satisfying the limiting absorption principle (|34[> and solutions v G H 2 (A4 R ) 
to the problem (|38| with F — fx and 3A ^ 0. In the next theorem we prove that this difference (error) 
decays with an exponential rate as R — > +oo. In other words, we prove exponential convergence of the PML 
method. The problem with finite PMLs can be solved numerically with the help of finite element solvers; 
discretization produces yet another error that we do not estimate here. 

Theorem 3. Let A G O a \ M. and let (3 be in the interval in (|35l) . Assume that /Jo 6 K \ er(A^f) is not an 
eigenvalue of the selfadjoint Neumann Laplacian °A in L 2 (A4 ) . Then there exists Rq > such that for any 
f G A and for all R > Rq a unique solution v R G Hq(A4 r ) of the problem (|38[) with F — fx converges on 
Ai r to the outgoing solution u_ (resp. the incoming solution u+ ) of the equation (°A — po)u — f in the case 
SA > (resp. 3 A < 0) in the sense that as R — > +oo the estimate 

(51) \\u ± -v R \\ 2 H2{Mr) ^Ce- 2 ^\\e^fx\\ 2 

holds with a constant C independent of R > Rq and f. 

Let us recall that the set A of analytic functions is dense in L 2 (A4). In particular, our construction and 
Theorem [3] remain valid for any / G L 2 {Ai) supported in Ai r - 
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Proof. By Theorem [T] it suffices to prove the estimate (jS"Tj) with u± replaced by u\. Let p E C°°(R) be a 
cutoff function such that p(x) = 1 for x ^ and x( a: ) — f° r x ^ 1/2- We set pr{x) = p(x — R) on II and 
pji = 1 onM\II. Then prU\ — Vr — u\ — vr on A^_r and the difference pru\ — vr E H 2 {M. r) satisfies 
with F = (,0^ — 1)/a + [ A A,/9_r]ua- Observe that 

\\(QR-l)fx\\ L * { Mn)<Ce-e R \\e( 3s h\\, \\{ x A,pR}u x \\ LHMR) ^Ce-P R \\e^u x \\ HHM) . 

This together with (f55"j) gives 

(52) ||F|| i2(A4j?) <ce-^||e^/A||. 
By Theorem [2] we have 

||pijUA - v R \\ H 2 iMR) < C||i ;l || L 2(_ Mi?) , i? > i? . 

This together with 

||«± ~ Uflll i» (Mr) ^ \\PR U \ - v r\\hz(Mr), R> Ro> r, 
and (|52l) completes the proof of (fS"Tj) . □ 
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